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Abstract 

We consider open systems where cars move according to the deterministic 



Nagel-Schreckenberg rules [24| and with maximum velocity Vmax > 1) what is 
an extension of the Asymmetric Exclusion Process (ASEP). It turns out that 
the behaviour of the system is dominated by two features: a) the competi- 
tion between the left and the right boundary b) the development of so-called 
"buffers" due to the hindrance an injected car feels from the front car at the 
beginning of the system. As a consequence, there is a first-order phase tran- 
sition between the free flow and the congested phase accompagnied by the 
collapse of the buffers and the phase diagram essentially differs from that for 

V,raax = 1 (ASEP). 

I. INTRODUCTION 

Driven diffusive processes have been widely studied as prototypes of non-equilibrium systems 
( 0^ 0)- Tliey are modelled as a lattice gas and are cliaracterized by a constant external 
force (e. g. electrical field) which sets up a steady current transporting information from 
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the boundaries to the bulk of the system. 

A well-known modification of the basic one-dimensional diffusive system is the asymmetric 
exclusion process (ASEP) which was first solved by by Derrida et al for open boundary 
conditions. The ASEP is defined as follows: Consider a one-dimensional lattice of L sites. 
Each site i (1 < i < L) is either occupied by a particle (tj = 1) or empty (tj = 0). A particle 
on site i has the probability p of hopping one site to the right if site i+1 is empty. At the 
left boundary of the system a particle is injected with probability a if i = 1 is empty. At 
the right boundary a particle on i = L is removed with probability /3. The ASEP can be 
divided into four classes according to the order in which to perform hopping, injection and 
removal: 

(a) random - sequential update( [5] — 0) 

(b) ordered - sequential update( 0, [p!^ ) 



(c) sublattice - parallel update( — [|T2|) 

(d) parallel update( - [IH) 

A detailed overview over all update types is given in P3 . 



An interesting feature of the ASEP is that phase transitions occur as a function of the model 
parameters. Usually there is a low density/high current phase and a high density/low cur- 
rent phase reminiscent to the "free flow" and the "jamming" states in vehicular traffic ||2CI|| - 
p^ . Being a cellular automaton the ASEP and its generalizations are well-suited to serve 



as simple models for traffic problems since efficient analytical and numerical techniques have 
been developed for their study. 

As it is common for traffic simulations we will use parallel update in the following because 
this is the most effective among the four update types and shows the best congruence with 
real traffic data In Fig la we reproduced the main results for the ASEP with parallel 



update: Based on investigations on global density, current, density profiles and correlation 
functions it turns out that there are two regimes, free flow and jamming, which are separated 
by the a=(3-lme {a: injection rate, f3: extinction rate). All parameters have in common that 



they do not depend on the extinction rate /3 (injection rate a) in the free flow (jamming) 
regime. 

Comparing the ASEP with real traffic, however, it is obvious that phenomena like acceler- 
ation and slowing down are not included in the model. Here, cars either do not move at 
all or move one site per time step. It can therefore be said, that they move with maximum 
velocity Vmax = 1. In order to get more reahstic results, Nagel and Schreckenberg introduced 
a model where cars are able to drive with different discrete integer velocities v, < v 



Another interesting feature of the parallel update procedure is that it induces additional 
short-range correlations compared to other updating procedures. An essential part of this 
paper will be therefore devoted to the investigation of short-range correlation functions (Sec- 
tion V) which have been already studied in corresponding systems with periodic boundary 
conditions for Vmax > 1 p6[ - and in systems with open boundary conditions for Vmax 



= 1 [0. Moreover, it turned out that correlation functions are well suited to describe the 
free flow - jamming transition - [^. 

The most significant difference between systems with open and periodic boundary conditions 
is the car density p. In a periodic system the car density is a tuning parameter and the 
probability to find a car at a certain site i is p. In systems with open boundary conditions, 
however, the situation is different as we have to deal with two different tuning parameters, 
namely the injection rate a and the extinction rate P and the density is a derived parameter. 
The influence of a and (3 on the car density implies that quantities like global density, cur- 
rent, and the density profile, which were studied for the ASEP {vmax = 1) in Il25|, ||31[ - 



|3^ , show a different behaviour than in periodic systems. For the case Vmax > 1 and open 
boundary conditions, however, only few results exist. Therefore, the cases Vmax = 1 and 
Vrraaz > 1 in systems with open boundary conditions will be compared with each other in 
this paper, too. 

The paper is organized as follows: In the next section the model is described. The current 
and the global density of the system are considered in Section III, in particular for the cases 



(3 = 1, a = 1, and l3 = 1-a. In Section IV we analyze the corresponding density profiles 
and in Section V the short-range correlation functions. The results are summarized and 
discussed in Section VI. 

II. MODEL 

Our investigations are based on a one-dimensional probabilistic cellular automaton model 
introduced by Nagel and Schreckenberg ||2^ . According to the Nagel - Schreckenberg (NS) 
model, the road is divided into L cells of equal size and the time is also discrete. Each site 
can be either empty or occupied by a car with velocity v = 0, 1, . . ., Vmax- All sites are 
simultaneously updated according to four successive steps: 

1. Acceleration: increase v by 1 if v < Vmax- 

2. Slowing down: decrease v to v = d if necessary (d: number of empty cells in front of the car). 

3. Randomization: decrease v by 1 with randomization probability p if p > 0. 

4. Movement: move car v sites forward. 

It is obvious that the NS - model is identical with the ASEP model with parallel update for 
maximum velocity Vmax = 1- In this paper, the randomization probability is p = 0, i. e. 
step 3 (randomization) is ignored. The investigations are mainly focused on Vmax = 5 but 
for comparison also maximum velocities Vmax = 2, 3, 4, 6, 7, ... are considered (see Section 
III). 

Open boundary conditions are defined in the following way: 

The system consists of L sites i with 1 < i < L (for the numerical simulations: L = 1024). 
At site i = 0, that means out of the system a vehicle with the probability a and with the 
velocity v = Vmax is created. This car immediately moves according to the NS rules. If i = 
1 is occupied by another car so that the velocity of the injected car on i = is v = then 
the injected car is deleted. At i = L+1 a "block" occurs with probability 1-/5 and causes a 
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slowing down of the cars at the end of the system. Otherwise, with probabihty /?, the cars 
simply move out of the system. 

III. CURRENT AND GLOBAL DENSITY 

The phase diagrams for systems with maximum velocities ^max = 2, 3, 5 are shown in Figs 
Ib-d. Fig lb resembles the case \max = 1 except for some deviations which are due to the 
fact that in systems with Vmax = 2 we do not have a particle-hole symmetry as for Vmax — 
1. The course of the free flow - jamming border for the case Ymax = 3, on the other hand, 
is very different (Fig Ic). Here, the a = /3 - line does not separate the free flow and the 
jamming regime. Instead, the jamming regime is larger than the free flow regime, and for 
high extinction rates P cars freely move for all a. For the maximum velocity Vmax — 5 these 
features are even stronger developed as it is obvious from Fig Id. 

Let us have a closer look at the /3 = 1 - line. The current q in Fig 2a increases with 
increasing a; for Ymax > 5 we have q{a < 0.5, /3 = 1) = a. For high injection rates, however, 
the curves surprisingly decrease (if Vmax ^ 4). This phenomenon cannot be observed in 
systems with maximum velocities Ymax = 2, 3 and for Vmax = 4 it is extremely weak. The 
maximum of the current is at a ~ 0.9 for Vmax = 5 and at a ~ 0.835 for higher maximum 
velocities. 

The corresponding global density p results from the current in Fig 2a by the relation 

p{a, p^l)^ 

^max 

as all cars freely move with maximum velocity Vmax- 

Considering the current (Fig 2b) and the global density (Fig 2c) for the injection rate a 
= 1, we see that for Vrnax — 2 these quantities behave similarly to the case Vmax — 1- For 
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'Vmax > 3 astonishing effects are observed which do not depend on the maximum velocity if 
^max ^ 5: Coming from low extinction rates /3 the current for Ymax ^ 5 increases propor- 
tionally to /3 and abruptly becomes constant at /3c = 0.835. For the global density, on the 
other hand, the transition seems to be continuous. 

Investigations of systems for large system sizes, however, show that the continuous change 
in the global density is just a finite size effect: Although the curves are qualitatively the 
same as those in Fig 2c the transition from free flow to jamming becomes more and more 
abrupt with increasing system size L. Furthermore, it turns out that the value of /3c is slightly 
smaller than for L = 1024. As a consequence from numerical investigations of systems with 
large L it is fair to assume that for L — >^ oo the current is described by 

5 4 
q{a = < -,ymax > 5) = -/3 jamming 

5 2 
q{a = 1, /3 > -, Vrnax > 5) = - free flow 

and the corresponding global density is given by 

_ 5 4 

p(q; = < -,v„ax > 5) = 1~5^ jamming 

_ 5 2 

p{a = !,/?> -,ymax > 5) = free flow 

For increasing system sizes current and global density converge to these values which can be 
"calculated" analytically as it is demonstrated in the following. Unfortunately, there exist 
no extensive analytical theory of the NS model for maximum velocities Vmax > 1- We must 
therefore restrict ourselves to a kind of bookkeeping which is nevertheless well-suited for the 
understanding of what is going on in the system. Furthermore, it should be emphasized that 
the representations of the conflgurations are snapshots between the slowing down step and 
the movement step. This is just a convention and does not change anything in the physical 
meaning. 



In order to get a better insight in the behaviour of the current and the global density 
we consider the special case a — P — 1. The car velocity is represented by numbers in 
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brackets, (v) = (0), (1), {vmax), and k connected unoccupied sites by the symbol x . 
The first number in brackets represents the car at i = where cars are injected. Then we 
have for 

t = : (Vmax) 

t = 1 : {ymax - 1) x^^'^^-^ (vmax)a:^'^~^'""" 

^ — ^ ■ \^max ^) y^max) y^maxj 



Now, a pair of consecutive cars is focussed found at the beginning of the system at time t 
= n with 2 > n > Vmax - ^■ 

t = n : {vmax - n) x^---" (v^„^ - n + 2) ■ ■ • (vmax) x^-'^^-^'^^. 

The difference of the velocities of the cars is Av = \ front - Vfeacfc — 2 and the velocity of each 
car increases by 1 due to the acceleration step of the NS model. Consequently, the space 
between the cars grows with Av t = 2t. After n — 1 time steps we have 

t = 2n - 1 : ■ • ■ (v„„, - 1) (v^„,) • • ■ (v^„,) ^^^-(^^-Dv^.. 

and finally, we find 

+ — On ■ l\T \ ^^max+Tl-l ^ . . . (it ^ max 

^ — ■ y^max) y^max) \^ max ) 

From now on, the space between the front and the back cars keep constant and consists of 
maximally 2(v^aa; — 1) empty sites due to n < Ymax ~ 1- 

The situation is different for the case n = Ymax- 

t = ymax- (0)(2)---(v„„,)x^-^— 

According to the left boundary conditions the car at site i = with velocity v = is deleted 
and a new car is created instead at the next time step: 



t = Vma. + 1 : (2)x2(3) • • ■ (v^„,)x^-v— 

Here, Av — 1 and the space between the cars grows as Av t = t. After Ymax-'^ time steps 
we finally get 

+ — O-fr 1 . {-IT \rf^^max f-tT \ . . . { -^r N ,>^^~Vmax (2VTnax — 1) 

^ — ^''max ^ • \^max)-^ \^ max ) y^'maxj-^ 

If one proceeds, it can be clearly seen that there are three scenarios (m = 0, 1, 2, • ■ •): The 
car created at site i = and t = n (a) is deleted according to the left boundary conditions 
if n = Vmax + 3m. (b) has Vmax empty sites in front if n = Vmax + 1 + 3m (c) has 2(yrnax — 1) 
empty sites in front if n = Vr„aa; + 2 + 3m. In other words, a self- repeating pattern establishes 
itself after a while according to 

{Q\ ... i\T \ rJ^^ fnax — ^ {-.J \ max { -.j \ ™2(Vmaa:~l) { -.j \ max {-.j \ r^^^ max 

y'max) y^maxj \^max) \^max) •'^ y^maxj 

/' 1 ^ rf^- ... f-iT \ rp'^O^max 1) ( -.T \ max ( -xr \ rY^'^^y max 1) ( -tr \ max f-fr \ ^^(^max 

{,-'-) \^max) •'^ \^max) \^ max j \^ max j \^maxj 

('}\ t"^ ■ ■ ■ (\r \ ^2(Vmax-l) (^r \ r^max (^j \ ^2(Vmaa:-l) (^j \ ryS max (^j \ ^2(Vm 

\y) yi^) jy y^maxj y^maxj y^maxj y^maxj y^maxj 

This is perhaps astonishing because we naively would expect Ymax unoccupied sites between 
two neighbouring cars for a = 1. Actually, there are also spaces consisting of 2{Ymax-^) 
sites what is a consequence of the hindrance the injected cars feels from the front car at 
the beginning of the system. In other words, Vmoa;-2 additional sites - so-called "buffers" 
(the motivation for this name will be explained later) - occur playing an important role for 
systems with maximum velocity Ymax > 3 as we will see below. 

Besides, our reflections clearly show that one has to wait for at least t = Y^ax time steps 
until the self-repeating pattern is established. Within this time period the first car created 
at t = has moved onto site i = v^^^,. Therefore, our considerations are only valid for 
systems which size is much larger than v^^^,, otherwise, boundary effects must be taken into 
account. 

From the self-repeating pattern follows that the distance between two neighbouring cars 
driving with Y^ax is alternately di = Ymax and d2 = 2(v^aa; - 1), i.e.. 



di 


= 2, d2 


= 2 


for 


^max 


= 2 


di 


= 3, d2 


= 4 


for 


^max 


= 3 


di 


= 4, d2 


= 6 


for 


^max 


= 4 


di 


= 5, d2 


= 8 


for 


^max 


= 5 


di 


= 6, d2 


= 10 


for 


^max 


= 6 



That means, buffers occur only for Vmax ^ 3. 

^max = 2 is a special case behaving similarly to Vmax = 1- It is therefore no surprise that the 
corresponding phase diagram, the global density, and the current resembles the case Vmax — 
1. If finite size effects are left out of consideration the current is obviously given by 

2 

q{a = f3 = l,ymax > 1) = - 

and the global density by 

2 

■p(a^ l,Vmax > 1) = ^ 

what coincides with numerical results. 



We will now investigate the effect of the buffers for the extinction rate P — 1. For that 
purpose we consider a slightly smaller injection rate by working a "disturbance" in the a 
— (3=1 pattern, i. e., at each time step except for one a car is created at i = 0. As the 
self-repeating pattern consists of three time steps we have three possibilities to place the 
disturbance. In Appendix A the effect is illustrated for systems with Vmax > 5 (because 
of lack of space we set v = Vmax in Appendices A, B). It turns out that the movement of 
the cars does not change at all for possibility (c). For (a) and (b), however, the distur- 
bance influences the system for three time steps as three cars show a deviating behaviour in 
Appendix A. Having a closer look on the sites affected by the disturbance we see that the 
current qdi^tiP = 1, Vmax > 5) | and the global density PdistiP = 1, "^max > 5) = are 
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higher there than ior a — P — 1. As altogether Avmax sites are concerned by the disturbance 
the effect increases with increasing Ymax- 

Considering the site i = in Appendix A it is obvious that the effect of the disturbance is 
different for maximum velocities Vmax < 5 as cars driving with Vmax — 4 cannot be injected 
with V = 5, cars driving with Vmax — 3 not with v 4 and so on. We do not go into details 
but just list up the results: Placing a disturbance at the beginning of a system with Ymax — 
2, 3, 4 one gets 

12 3 

(a) (IdistiP = 1, "^max = 4) = — ' PdistiP = i,ymax = 4) = — 
^distiP = 1> "^max = 3) = ^, PdistiP = 1> "^max = 3) = ^ 

2 _ 1 

^distiP = 1> "^max = 2) = -, PdistiP = 1, V^ac, = 2) = - 

3 3 

(b) (IdistiP = 1> "^max = 4) = ^, PdistiP = 1, "^max = ^) = 

no effect of disturbance for Vmax — 3 

^distiP = 1> '^max = 2) = ^, PdistiP = 1, ^max = 2) = | 

(c) no effect of disturbance for Vmax = 2, 3, 4 

Superposition of all possibihties (a), (b), and (c) leads to the result that the effect of the 
disturbance is weaker for Ymax = 4 than for corresponding systems with Vmax > 5. For 
maximum velocities Vmax = 2, 3 the current and the global density decrease and that is why 
the maximum of the curves in Fig 2a is at a = 1 if Vmax < 3. 

As far as the position of the maximum of the current is concerned we can only give a 
hand- waving argument: It is obvious from Appendix A that the disturbance affects the 
development of two buffers. On the other hand, it can be easily seen that ior a — f3 — 1 a 
buffer is created every three time steps (and consequently, two buffers are created in six time 
steps) . Therefore, the strongest effect is expected when the system is disturbed with the 
rate (l-cc) = |. If (l-a) becomes higher the buffers being necessary for the increase in the 
current and the global density cannot develop. This may be the reason why the maximum 
for the curves in Fig 2a with Ymax > 5 is found at ct ~ |. 
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For the injection rate a = 1 the buffers have an even more dramatic effect which can 
be observed at the end of the system. In analogy to /3 = 1 we start with the special case a 

— (5 — 1. By simple analytic considerations it turns out that a self-repeating pattern 

. , . ^2(VTnax 1) (-IT \ max (-it \ /y.2(Vmax~l) (-it \ max (-it \ rj^^ max 1) (-tj \ 

\^ max j •'^ \^ max ) \^ max j \^ max ) \^ max ) 

{ \r \ rf^^ max 1) ( -tr \ max ( -tr \ max 1) (-tr \ max (-tr \ max 1 

V * max ) \^ max J \^ max / V max / **' V max ) 

(-tr \ max ('IT \ rp^i^max 1) f \ rp^max ('it \ rp^O^max ^^{\r \ rp^ max 

^ max ) \^ max ) \^max) \^ max ) \^max) 

establishes itself at the end of the system, too (L v^^^) . It is important to mention that 

- due to (3 — 1 - no blockage occurs at all at the right boundary and that the buffers reach 
the right boundary with the rate abuffer = 5- The introduction of a disturbance (i.e. the 
consideration of an extinction rate which is slightly smaller than /3 = 1) means here to place 
a single blockage at the end of the system. According to the self- repeating pattern consisting 
of three time steps we have to consider three possibilities. Prom Appendix B it turns out 
that the ^max-'^ additional sites (resulting from the hindrance the cars feel at the beginning 
of the system from the front car) play an important role at the end of the system, too. Here, 
they have the effect of a "buffer" against the influence of the right boundary. It can be seen 
from Appendix B that two buffers are necessary to neutralize the blockage effect at the end 
of the system. Therefore, as long as (1-/3) < ^Oibuffer — ^ ^ jamming wave cannot develop. 
For /3c = |, however, there is a jump in the global density (remember that our analytical 
considerations are based on systems with size L — ^ 00, for L = 1024 the change from free 
flow to traffic is less abrupt due to finite size effects). As mentioned above we have p{a — 
1, P > f3c, Vmax > 5) = 3^p^ in the free flow regime and 'p{a = 1, /3 < /3c, Vmax > 5) = 
1 - 0.8/3 in the jamming regime. At /3c = | p immediately increases from — (free flow) 
to I (jamming). That means that there is a jump of '^y"^"^ at the critical extinction rate 
what corresponds to the buffer density in the free flow regime. In other words: At /3c — 
I the buffers cannot neutralize the blockage at the right boundary any longer. The buffer 
effect breaks down, jamming waves propagate from the end of the system to the left, and the 
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buffers (vmax - 2 sites on Svmax sites each) are completely occupied by cars. Consequently, 
both current and global density show similar behaviour as the corresponding quantities for 

Vrnax = 1 if < |. 

Another interesting feature observed in Figs 2b, c is that current and global density do 
not depend on the right (left) boundary conditions, i.e. not on /3 (not on a and Vmax), if 
the system is in the free flow (jamming) regime. This is not only valid for a = 1 but also 
for general injection and extinction rates as it can be seen in Fig 3a for the current and for 
Fig 3b for the global density. 



In order to compare our results with those for corresponding systems with periodic boundary 
conditions we investigate the case f3 = 1-a. For f3 = 1-a there are rather similar conditions 
at the left and at the right boundary and therefore, systems with open and with periodic 
boundary conditions can be compared at best with each other. 

The fundamental diagram for systems with periodic boundary conditions (PBC) is com- 
pletely determined by the maximum velocity Vmax (see [^| and references therein). The 



current of the system is given by q^^'~^{p < pc) = VmaxP for freely moving and by q^-^'"(p > 
Pc) = I - p for jammed cars. The critical density is given by p^ = y ^ 
In the case of open boundary conditions, on the other hand, it turns out from numerical 
results for Vmax > 5 that the current in the free flow (jamming) regime only depends on the 
injection (extinction) rate according to 

q{f3 = 1-a) = a for a < ac, (3 > Pc 
q(/3 = 1-a) = 0.8(3 for a > ac, (3 < Pc 

and consequently, the transition takes place at Oc = 0.44 {(3c = 0.56). The global density for 
(3 = 1-a shows finite size effects as in the case a = 1. For large L, however, the transition 
from free flow to jamming becomes sharp. Then the global density is described by 
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V, 



max 



a 



for a < ac, P > Pc 



p{P = 1-a) = 1 - 0.8/? 



for a > ac, P < Pc 



with a jump at a. 



■c 



4 
9 



0.44 {P, 



5 
9 



0.56) 



The results for P — 1-a induce the identity 



q(/3 = 1-a) = q(/? = 1) 



for a < ac, P > Pc 



q{P — 1-a) — q{a — 1) 



for a > ac, P < Pc 



This indicates that the movement of the vehicles in the high density or jamming regime is 
dominated by the right boundary conditions, in the low density or free flow regime by the 
left boundary conditions. To get a better insight in this question we will have a closer look 
on the density profiles and the short-range correlation functions which are analyzed due to 
the three special cases (see also Fig Id) 

1. P = 1: shows the influence of the left boundary 

2. a — 1: shows the influence of the right boundary 

3. P = 1-a: systems with open and periodic boundary conditions can be compared at best 
with each other 

The investigations of this section clearly show that the case Vmax = 5 includes all features 
which are characteristic for higher maximum velocities, too. For this reason we conflne our- 
selves to systems with Vmax — 5 (and L = 1024) in the following. 



IV. DENSITY PROFILES 
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A. p = 1 



In this section we investigate the influence of the left boundary on the density profiles. The 
best way to do this is to set /? = 1, because in that case the right boundary has no influence 
on the system. 

From Figs 4a,b it can be seen that the density profiles are characterized by a periodical 
structure. This is a significant difference to the case Vmax = 1 where oscillations cannot be 
found at all For Vmax = 5, however, the density profiles show a certain pattern recurring 
with the period Ai = 5. In order to understand this phenomenon we consider the density 
profiles for very low injection rates first. 

For a = 0.05 (see Figs 4a,b) the probability of generating a car at i = in two successive 
time steps is very small and therefore, the cars at the beginning of the system do not interact 
with each other. That means that a car which is created on i = with velocity 5 (according 
to the left boundary conditions) moves to i = 5 at the next time step and can be found on 
the site i = 5n after n time steps (n = 1, 2, 3, ...). The density on these sites is p ~ a. As 
it is obvious from Figs 4a,b a car can be also found on i = 5n+4 for small a, too, but the 
probability for that is very small. 

For increasing injection rates a, however, the probability of generating cars in two successive 
time steps increases and with it the hindrance a car at the beginning of the system feels 
from the front car. This can be understood as follows: Let us create a car A at time step t 
and a car B at time step t+1. Considering the system at t+1 we see that car A is on i = 5 
having the velocity 5 whereas car B on i = has the velocity 4 because there are only four 
empty sites to car A. At the time step t+n, car A is on i = 5n and car B on i = 5(n-l)-l. 
To sum it up it can be said that the hindrance due to the left boundary conditions leads to 
a shift of the position of the cars within the system. This shift is reflected in the periodic 
pattern of Figs 4a,b. Whereas it is rather probable to find a car on i = 5n+5 and on i = 
5n+4, the probability of finding a car on i = 5n+2 is much smaller and for i = 5n+3 it 
nearly vanishes. 
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The most important result, however, is the fact that the sites i = 6+5n are never occupied 
according to the left boundary conditions so that the density on these sites have the value 
p(i=6+5n) = for all a. Before turning back to this point we have a look at the case a ~ 
(3 — 1 which is of special interest in the following section, too. 
For a — (3 — 1 the corresponding density profile has the following form: 

p{\) = 1 if i = 5n+4 or i = 5n+5 
p(i) = else 

as it can be easily deduced from the left boundary conditions. 

B. a = 1 

We investigate the influence of the right boundary now. Unfortunately, the influence of the 
left boundary cannot be completely left out of consideration by setting, for example, o; = 0, 
because in that case no cars would be generated at all. Instead, we choose a = 1, because 
only for q; = 1, the cars are deterministically created. This allows us to distinguish between 
the influence of the right and of the left boundary. 

In Fig 5a we can sec that the situation for a = 1 is very different from that described in the 
previous section. For high extinction rates we still recognize the periodic structure already 
known from the case (3 — 1. For extinction rates (3 between 0.75 and 0.85 something inter- 
esting happens: the oscillations vanish and the envelope of the density profile rises. For low 
extinction rates the density profiles are just a constant which value increases with decreasing 
(3. 

In order to understand this change we consider density profiles for 0.83 < (3 < 0.84 in Figs 
5b-d. On i = 4-|-5n and i = 5+5n we find p(i) = | resulting from the infiuence of the left 
boundary (see Section IV.A). The other sites (with p — lor (3 — 1), however, increasingly 
reflect the influence of the right boundary with decreasing extinction rates. Coming from 
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high (3 the density on i = 6+5n, i = 7+5n and i = 8+5n seems to "come away" from the 
p(i) = - hne starting at the right boundary. This phenomenon can be explained due to 
the repulsion the car feels at the right boundary with decreasing probability (3 of being ex- 
tincted. Consequently, a jam develops at the right boundary which expands to the left. For 
(3 ^ 0.837 the influence of the right boundary finally reaches the beginning of the system 
(Fig 5c). For (3 = 0.84 the sites i = 4+5n and i = 5+5n indicate the repulsion at the right 
boundary, too, as the density profile becomes p > | there. In parallel to this the oscillations 
resulting from the left boundary conditions vanish, a process which starts from the end of 
the system as well. 

Our observations have been quite qualitative so far. In the following the transition de- 
scribed above will be analyzed in detail and for that purpose we will have a closer look at 
the sites i = 6-|-5n. As we know from the previous section these sites are never occupied 
according to the left boundary conditions. In other words: The occupation of the sites i — 
6+5n is exclusively an effect of the right boundary. Therefore, these sites play an important 
role as they show the repercussion of the right boundary on the system. 
The density on these sites is shown in Figs 6a-c. The density profiles correspond to the same 
(3 as in Fig 5b but here, all sites except for i = 6+5n are left out of consideration. Let us 
first consider the density profiles for (3 > f3c (Fig 6b) which are exponential functions p (i 
= 6-|-5n) = Pmax{ P ) e*^^^^(^"^) {pmax{ )'■ maximum value of the density on the sites i = 
6+5n). In Fig 6d the exponent c(/3) is drawn against the extinction rate (3 and it is obvious 
that c(/3) — 'k {(3 - (3c) with /3c = 0.8362 and k Ri 2. Whereas /3c can be clearly identified as 
the critical extinction rate where the transition from freely moving to jammed traffic takes 
place the factor k is still an open question. 

If we pass over to the density profiles for (3 < (3cii can be easily seen in Figs 6c,d that the 
density profiles have the form p(i = 6-|-5n) = Pmax (/?)[! - 

The behaviour of the density profiles described in this section has the following physical 
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explanation: 

As it is well-known, the right boundary has no effect on the density profiles for /3 = 1. With 
decreasing /3, however, there is a growing probability of a blockage at the end of the system, 
i. e., cars are increasingly hindered from moving out of the system. Consequently, a jam 
develops showing the growing influence of the right boundary with decreasing /3. For /5 > /3c 
the influence of the right boundary exponentially diminishes (Fig 6b). Fig 6a further shows 
that the left boundary conditions are still valid for the whole system, what can be seen at 
the oscillations of the density profile and in the constant value p (i) = | on the sites i — 
4-|-5n and i = 5-|-5n characteristic for the case a = P = 1. For decreasing /3 the jam and 
with it the influence of the right boundary expands to the left. 

At /3 — /3c the repercussion of the right boundary reaches the left boundary, and the decay 
of the jam is proportional to i. Simultaneously, the influence of the left boundary is still 
present in the whole system, too, which manifests itself in the oscillations in the density 
profile going from the left to the right boundary (Fig 6a). So it can be said that for the 
extinction rate (3 = the influence of the left and that of the right boundary coexist in the 
whole system. 

However, beginning from the right the oscillations vanish when the extinction rate is further 
decreased (Fig 6a). This indicates that the influence of the left boundary is pushed back 
for P < (3c- The form p(i = 6-|-5n) = Pmax(P)[i - e^^*^) ^] shows the decrease of unoccupied 
sites and may be a hint at a symmetry around the transition point. 

For very small (3, the left boundary does not have any relevance at all for the movement of 
the cars in the bulk. 

Finally, let us say some words about the maximum value Pmax{/3) = max[p{i = 6+5n)]. 
From Fig 6a it is obvious that Pmax{P) can be identified with the density on site i = 1021, 
Pmax{(3) — p {P, i — 1021). From Fig 6e it turns out then that 

PmaxiP) = PmaxiPc) e^'^''-''^^ for /3 > 
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(ki -24.46; k2 ^ -0.8; /5c = 0.8362). Thus, the transition from freely moving to jammed 
traffic is reflected at the right boundary, too. 

C. /3 = 1-a 

We have already mentioned that for P — 1-a we have rather similar conditions at the left 
and at the right boundary and therefore, systems with open and with periodic boundary 
conditions can be compared at best with each other in that case. 

We must keep in mind, however, that there are significant differences for /3 — 1-a, too, 
especially if the randomization probability is p = 0: In systems with periodic boundary con- 
ditions the movement of the cars is fully deterministic and the car density p in the system 
keeps constant. Each site in the system has the same probability of being occupied and 
therefore, the density profiles of systems with periodic boundary conditions are constants 
with the value p (the latter statement is also valid for randomization probabilities p > 0). 
For corresponding systems with open boundary conditions - due to the injection rate a and 
the extinction rate 13 - we always have a non-deterministic element at the boundaries of the 
system, also for the randomization probability p = (which only refers to the movement in 
the bulk). 

Generally speaking, the density profiles for /3 = 1 - a show a similar behaviour as those 
for the case a — 1: For very low extinction rates (and high injection rates) the density pro- 
files are identical with the density profile of a corresponding system with periodic boundary 
conditions. For high /3 (and low a) the density profiles show the periodic structure already 
known from the previous sections as a typical feature of the free flow regime. At (5c — 0.56 
(and ac = 0.44) the transition from free flow to jamming takes place. For (3 > (3c the curves 
have the form p (i = 6+5n) = p^„,( (3 ) eC(/3)(i-L), for p < Pc p(i = 6+5n) = PmaxW - 
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and for /3 = /3c we have a straight hne. The only difference to the a = 1 case is the 
value of the critical extinction rate and of k: For (3 = 1-a we have pc = 0.56 and k = 3.75. 

In Section III we have already mentioned that in the high density regime the global density 
(current) for = 1-a is identical with the global density (current) for a = 1 and in the low 
density regime with the global density (current) for (3 = 1. From Figs 8a,b it is obvious 
that similar effects can be also observed for the density profiles, too. Having a closer look 
at them we see that the profiles for (3 = 1-a and (3 = 1 are identical if the injection rate a is 
low. For increasing a the density profiles for (3 = 1-a start to lift at the end of the system 
indicating the growing influence of the right boundary on the system for increasing a. On 
the other hand, comparing the density profiles for (3 = 1-a and a = 1 with each other we 
see that they are identical for very low (3. For increasing P the density profile "drops" at 
the beginning of the system. Accordingly, this behaviour shows the growing influence of the 
left boundary on the system. In the transition regime, however, the density profiles for (3 = 
1-a are very different from those for the cases a = 1 and (3=1. 



V. CORRELATION FUNCTIONS 

In this section we consider correlation functions 

C(i,t) = < r/(i',t'), 77(i'+i,t'+t)>i.^t' -P' 

for short ranges with 

ri{i\V) = 1 if site i' is occupied at time t' 
r/(i',t') = else 

This kind of correlation functions has been already investigated for systems with periodic 



boundary conditions and the randomization probability p = 0.5 in It turned out that in 
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the free flow regime tliere is a propagating peak at i = Vmaxt with a shoulder at i = Vmaxt - 1 
and with anticorrelations around it. The density where these anticorrelations are maximally 
developed is defined as the density where the transition from free flow to jamming takes 
place. For higher densities a jamming peak occurs at i = -1 



It would be interesting to see if these features can be also found for systems with open 
boundary conditions. But considering the deterministic case in this paper we should inves- 
tigate the correlation functions for systems with periodic boundaries and p = flrst. From 
Fig 8a it can be seen that the propagating peak is sharp and that there are further peaks at 
i = Ymax^ ± 6n (n = 1, 2,...) as the movement of the cars in the ring is deterministic. Due 
to the fact that the initial conflguration is random, however, these peaks become smaller 
and smaller with increasing n. Between the peaks anticorrelations are observed which are 
best developed around the peak at i = v^axt. Generally speaking it can be said that in the 
free flow regime the correlation functions C(i,t) are symmetric around the site i = 'Vmax^- 
Coming from low densities the anticorrelations become deeper and deeper with increasing p. 
At p = Pf. = ,r 1 1 the car distribution is well-deflned: all vehicles drive with the maximum 
velocity ^max = 5 and between two neighbouring cars there are ^max = 5 empty sites each. 
Correspondingly, the correlation function for p = pc is periodic with 

C(i,t) = p -p^ if i = ^maxt ± 6n 
C(i,t) = - p2 else 

At this density where the transition from free flow to jamming takes place the anticorrela- 
tions reach their minimum. 

For higher densities a jamming peak develops at i = -1 (due to the hindrance the back car 
feels in the jam) with anticorrelations at i = ± 1. At all other sites peaks and anticorrela- 
tions vanish. If the density further increases fewer and fewer cars move (with v > 0) and 
therefore, the anticorrelations at i = ± 1 disappear. Corresponding to the symmetry around 
i = Vmaxt in the free flow regime, the correlation functions for p > Pc are symmetric around 
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i = -1. 



Let us turn back to systems with open boundary conditions which is the real topic of this 
paper. In Fig 8b we consider correlation functions from the middle of the system because 
the influence of the boundaries is minimal there. It is obvious that for high densities the 
correlation functions in systems with open boundary conditions are nearly identical. Merely 
the minor maxima at i = ± 2 in Fig 8a shift onto i = ± 3 in Fig 8b. 

If the density in the system is low, however, the situation is completely different: For sys- 
tems with open boundaries we have a random element at the boundaries where cars are 
randomly created and deleted at each time step. Therefore, due to the permanent presence 
of randomization even if the movement in the bulk is deterministic we can only observe the 
propagating at i = v^axt (and a very small one at i = Ymax'^ =t 6). Around the propagating 
peak there are anticorrelations, too, but they are not so well-developed as the anticorre- 
lations of corresponding correlation function in the case of periodic boundary conditions. 
However, a common feature of systems with open and periodic boundary conditions is the 
symmetry of the correlation functions around i = ymax^ in the free flow regime. 

As we have already mentioned the anticorrelations around the propagating peak play an 
important role in systems with periodic boundary conditions and we will now discuss the 
question if similar features can be observed for systems with open boundary conditions. In 
Figs 8b-d we consider short-range correlation functions at the beginning, the middle, and the 
end of the system. Coming from high extinction rates (3 (with /3 — 1-a) the anticorrelations 
become deeper and deeper everywhere in the system. But what it is interesting is the fact 
that the anticorrelations start to vanish again at an extinction rate where the influence of 
the right boundary reaches the corresponding sites (see also Section V.C). Strictly speaking 
the anticorrelations start to vanish at about P — 0.42 (at the end), /? = 0.43 (in the middle), 
and f3 = 0.44 (at the beginning). 

Therefore to sum it up it can be said that the injection rate (extinction rate) where the 

21 



probability to find a car in the neighbourhood of another car is minimal (that means, where 
the anticorrelations start to vanish) can be considered as the injection rate etc (extinction 
rate (3c) where the transition from free flow to jamming takes place. 



VI. CONCLUSIONS AND DISCUSSION 

Systems with open boundaries where cars move deterministically with maximum velocity 
^max > 1 show interesting features mainly resulting from the competition of the left and 
of the right boundary for the influence in the system and from the existence of so-called 
"buffers". 

The latter plays a fundamental role at the comparison of systems with Ymax ^ 3 and Vmax 
— 1. One of the most important questions in this context is why the border between free 
flow and jamming for Vmax > 3 has such a different course than the corresponding border for 
the case Ymax = 1- By simple analytical considerations it turns out that - as a consequence 
of the hindrance an injected car feels from the front car - spaces > Ymax develop for high 
injection rates a (for the special case oi a — P — 1 there are alternately 2{\rnax-^) a-nd Vmax 
sites between neighbouring cars for all Ymax > !)■ That means in addition to the expected 
^max sites further sites occur which are the reason why the maximum current is found at 
a < 1 and | < /? < 1 for Ymax > 5. We call these additional sites "buffers" because they 
also have a buffer effect at the end of the system: Due to the buffers the development of 
jamming waves is suppressed up to an injection rate = | (for high a and Ymax > 5) and 
this buffer effect is responsible for the characteristic course of the free flow - jamming border 
for Ymax > 5. The transition from the free flow to the congested phase is of first order and 
accompagnied by the collapse of the buffers. 

In this context, it should be emphasized that the occurrence of buffers - and consequently 
the specific features of the Ymax > 3 model - is due to the parallel updating mechanism 
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and not an effect of the particular injection rule. Naturally, there are other possibilities of 
generalizing the ASEP to Ymax > 1, for example, one could keep the existence of the car at i 
= if i = 1 is occupied by another car. Simulations based on this alternative rule show that 
the phase diagram and the Q;,/?-dependence of the current are qualitatively the same as the 
corresponding quantities considered in this paper. This has been confirmed by analytical 
investigations of the special case Ymax — 5, a = /3 = 1 (according to Section III) where 
buffers occur, too. 

As global density and current (from now on we exclusively refer again to the injection rule 
defined in Section II) show no qualitative differences for Ymax > 5, a detailed analysis of 
the influence of the boundary conditions on the system (by means of density profiles and 
short-range correlation functions) is confined to the maximum velocity Vmax ^ 5. Further- 
more, our numerical investigations are based on systems with size L = 1024. It must be 
mentioned here that finite size effects occur which do not play an important role, however: 
the discontinuous transition in the global density becomes continuous and the value for the 
critical extinction rate {(3c = 0.836 for a = 1 and L = 1024) is slightly higher than in the 
case L — > oo. 

If we consider a single site of the system there are three possibilities: 

A. the site is exclusively under the influence of the left boundary 

B. the site is exclusively under the influence of the right boundary 

C. the site is under the influence of both the left and the right boundary 

In the free flow regime the system consists of sites of the types A and C, in the jamming 
regime of sites of the type B and C. The critical injection rate (the critical extinction rate 
Pc) where the transition from freely moving to jammed traffic takes place, is the only a {(3) 
where all sites of the system belong to the C-type, that means, where the influence of the 
left and of the right boundary coexists in the whole system. The farther we go away from 
the transition point the stronger is the dominance of the A-sites in the free flow regime and 
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of the B-sites in the jamming region. 

In the free flow (jamming) regime the current does not depend on the injection rate a (ex- 
tinction rate /3) what confirms the dominance of the left (right) boundary influence in the 
free flow (jamming) regime. 

Comparing the density profiles for Y^ax — 5 with those for \rnax = 1 we see that the most 
significant differences are found in the free flow regime: In the free flow regime the density 
profile for Y^ax — 5 shows periodic structure with the period Ai = Ymax — 5 which is due to 
the free movement of the cars. Another interesting result is the fact that the sites i = 6-|-5n 
(n = 1, 2, ...) are never occupied when they are beyond the sphere of influence of the right 
boundary. In the jamming regime, however, the density profiles for Ymax — 5 and ^max — 1 
are nearly the same. The investigation of the density profiles, especially the behaviour on 
the sites = 6 + 5n enables the precise localization of the phase transition. 

The short-range correlation functions C(i,t) show that there are parallels between systems 
with open and with periodic boundary conditions, which are the following: In the free flow 
regime C(i,t) is symmetric around the site i — v^axt and in the jamming region around i 
— -\ which may be a hint at a symmetry. Free flow (jamming) is characterized by a peak 
at i = v^aa;t (i = -1) with anticorrelations around it. Furthermore, in systems with open or 
periodic boundary conditions the anticorrelations around the free flow peak are maximally 
developed when the transition from free flow to jamming takes place. 
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VIII. APPENDIX A: DISTURBANCE AT THE BEGINNING OF THE SYSTEM 
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IX. APPENDIX B: DISTURBANCE AT THE END OF THE SYSTEM 
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FIGURE CAPTIONS 



Fig la: 

Phase diagram with density profiles for ^max = 1 in dependence on the injection rate a and 
the extinction rate /9 (according to [17]-[22]) 

Fig lb: 

Phase diagram in dependence on the injection rate a and the extinction rate /? {ymax — 2) 
Fig Ic: 

Phase diagram in dependence on the injection rate a and the extinction rate /9 {wmax — 3) 
Fig Id: 

Phase diagram in dependence on the injection rate a and the extinction rate /3 {ymax — 5). 
Our investigations are focused on the cases /3 = 1, a = 1, and /3 = l-ct marked by dashed 
hues. 

Fig 2a: 

Current for /3 = 1 and the maximum velocities \max = 2, 3, 10 
Fig 2b: 

Current for a = 1 and the maximum velocities ^max = 2, 3, 10 
Fig 2c: 

Global density for a = 1 and the maximum velocities \max — 2, 3, 10 
Fig 3a: 

Current in dependence on the injection rate a and the extinction rate (5 i^max = 5) 
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Fig 3b: 

Global density in dependence on the injection rate a and the extinction rate j3 {vmax — 5) 
Fig 4a: 

Density profiles at the beginning of the system (/5 = 1) 
Fig 4b: 

Density profiles at the end of the system {j3 = 1) 
Fig 5a: 

Density profiles for a = 1 
Fig 5b: 

Density profiles for a = 1 around the critical extinction rate /3c 
Fig 5c: 

Detail from Fig 5b at the beginning of the system [a = 1) 
Fig 5d: 

Detail from Fig 5b at the end of the system (a = 1) 
Fig 6a: 

Density profiles from Fig 5b taking only the sites i = 6+5n into account {a = 1) 
Fig 6b: 

Logarithmic plot of the density profiles for /3 > (3c taking only the sites i = 6+5n into account 
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Fig 6c: 

Logarithmic plot of p^ax - /o(i=6+5n) for a = 1 
Fig 6d: 

Gradient of the density profiles in Figs 6b, c depending on the extinction rate /3 {a = 1) 
Fig 6e: 

Maximum value of the density profiles p(i=6+5n) on i = 1021 (a = 1) 
Fig 7a: 

Comparison of the density profiles for (3 = 1-a and /3 = I 
Fig 7b: 

Comparison of the density profiles for (3 = 1-a and a = 1 
Fig 8a: 

Correlation functions for systems with periodic boundary conditions 
Fig 8b: 

Correlation functions in the middle of the system for (3 = 1-a 
Fig 8c: 

Correlation functions at the beginning of the system for (3 = 1-a 
Fig 8d: 

Correlation functions at the end of the system for /3 = 1-a 
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